Geometry-Driven Shift in the Tomonaga-Luttinger Exponent of Deformed Cylinders 
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We demonstrate the effects of geometric perturbation on the Tomonaga-Luttinger liquid (TLL) 
states in a long, thin, hollow cylinder whose radius varies periodically. The variation in the surface 
curvature inherent to the system gives rise to a significant increase in the power-law exponent of the 
single-particle density of states. The increase in the TLL exponent is caused by a curvature-induced 
potential that attracts low-energy electrons to region that has large curvature. 
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Studying the quantum mechanics of a particle confined 
to curved surfaces has been a problem for more than 
fifty years. The difficulty arises from operator-ordering 
ambiguities which permit multiple consistent quan- 
tizations for a curved system. The conventional method 
used to resolve the ambiguities is the confining-potential 
approach 0,0). In this approach, the motion of a particle 
on a curved surface (or, more generally, a curved space) 
is regarded as being confined by a strong force acting nor- 
mal to the surface. Because of the confinement, quantum 
excitation energies in the normal direction are raised far 
beyond those in the tangential direction. Hence, we can 
safely ignore the particle motion normal to the surface, 
which leads to an effective Hamiltonian for propagation 
along the curved surface with no ambiguity. 

It is well known that the effective Hamiltonian involves 
an effective scalar potential whose magnitude depends 
on the local surface curvature 0, 0, Gl IE]- As a re- 
sult, quantum particles confined to a thin, curved layer 
behave differently from those on a flat plane, even in 
the absence of any external field (except for the con- 
fining force). Such curvature effects have gained re- 
newed attention in the last decade, mainly because of 
the technological progress that has enabled the fabri- 
cation of low-dimensional nanostructures with complex 
geometry d 0, 1, ©, M, [U 03, El, ■ From the theoret- 
ical perspective, man y in trig uing phenomena pertinent to 
electronicstates [H, M, M, \MM, IMJM electron 
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diffusion [23j , and electron transport [2ij25U26l. l27j have 
been suggested. In particular, the correlation between 
surface curvature and spin-orbit interaction 1281. 1231 as 
well as with the external magnetic field [3(| l3ll. |32| has 
been recently considered as a fascinating subject. 

Most of the previous works focused on noninteract- 
ing electron systems, though few have focused on in- 
teracting electrons (33| an d their collective excitations. 
However, in a low-dimensional system, Coulombic inter- 
actions may drastically change the quantum nature of 
the system. Particularly noteworthy are one-dimensional 
systems, where the Fermi-liquid theory breaks down so 
that the system is in a Tomonaga-Luttinger liquid (TLL) 
state [Hj]. In a TLL state, many physical quantities ex- 



hibit a power-law dependence stemming from the absence 
of single-particle excitations near the Fermi energy; this 
situation naturally raises the question as to how geomet- 
ric perturbation affects the TLL behaviors of quasi one- 
dimensional curved systems. Peanut-sha ped Cqo poly- 
mers [35l |36| and M0S2 hollow nanotubes [37| are exem- 
plary materials to be considered for studying TLL behav- 
iors; they are thin, long, hollow tubules whose radius is 
periodically modulated along the tube axis. Hence, the 
periodic surface curvature intrinsic to the systems will 
produce sizeable effects on their TLL properties if they 
exhibit one-dimensional metallic properties. 

In this Letter, we examined the effects of geometric 
curvature on the TLL states in quantum hollow cylinders 
with a periodically varying radius and demonstrate that 
the presence of a curvature-induced potential can yield 
a significant increase in the power-law exponent a of the 
single-particle density of states n(uA near the Fermi en- 
ergy Ep; i.e., n(oj) cx \hu — Ep\ a [38|,|39j]. The geometric 
conditions required for the shift in a to be observable 
are within the realm of laboratory experiments, which 
implies that our predictions can be verified with existing 
materials. 

We first considered noninteracting spinless electrons 
confined to a general two-dimensional curved sur- 
face S embedded in a three-dimensional Euclidean 
space. A point p on S is represented by p = 
{x{u l ,u 2 ),y(u l ,u 2 ), z(u l , u 2 )), where (u l ,u 2 ) is a curvi- 
linear coordinate spanning the surface and (x, y, z) are 
the Cartesian coordinates in the embedding space. Us- 
ing the notation pi = dp/du l (i = 1,2), we introduced 
the following quantities: gij = pi ■ Pj, hij — Pij ■ n, 
n = (pi x Pj)/\\pi x pj\\, where n is the unit vector 
normal to the surface. Using the confining-potential ap- 
proach @, 0] , we obtained the Schrodinger equation for 
noninteracting electron systems on curved surfaces as fol- 
lows: 
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du 



(1) 



where g = det(py), <? lJ = [40], and m* is the effective 
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FIG. 1: (color online) Schematic illustration of a quantum 
hollow cylinder with periodic radius modulation. 



mass of electrons. The quantities K = (hixh 22 — h 2 2 )/g 
and H = (#11/122 + 322/111 - 2g\ 2 h\ 2 ) j (2g) are the so- 
called Gaussian curvature and mean curvature, respec- 
tively, both of which are functions of (it 1 , u 2 ). The term 
H 2 — JC in Eq. (pj is the effective scalar potential induced 
by surface curvature. 

We next focused on a hollow tube with a 
periodically varying radius represented by p = 
(r(z) cos r(z) sin #,z) (see Fig. [T]). The tube radius r{z) 
is periodically modulated in the axial z direction as 



, . Sr Sr ( 2ir 

r(z)=r -- + -cost -z 
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where the parameters r and Sr are introduced to ex- 
press the maximum and minimum of r(z) as ro and 
ro — Sr, respectively. Because of the rotational symme- 
try, the eigenfunctions of the system have the form of 
6) = e me ip n (z). Thus, the problem reduces to the 
one-dimensional Schrodinger equation 



^2 r 
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1 d r d 
rf dz f dz 



j+(H 2 -IC) ip n (z) = Eip n (z), 

" (3) 

where f(z) = Vl + r' 2 , K = -r" /(rf 2 ), and H = (f 2 - 
rr")/(2rf 3 ) with r' = dr/dz. 

Equation ^ is simplified by using a new variable 
£ = £(z) = Jq f(rj)dr], which corresponds to the line 
length along the curve on the surface with a fixed 9. 
Straightforward calculation yields (l9j 



d£_ 2 



U n (0 



2m*a 2 E 

h 2 
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with U n {0 = (n 2 - \)a 2 /r 2 - r"V/(4/ 6 ), where r, r", 
and / are regarded as the functions of £ using the in- 
verse relation z = z _1 (^). In order to derive Eq. (|4]), we 
introduced the length scale, a, and then multiplied both 
sides of Eq. © with a to make the units of U n and e 
dimensionless. Notice that by the definition of £(z), U n 
is periodic with a period A = £(A) depending on ro and 
Sr (as well as A). 

Figure [5] shows the spatial profile of U n within one 
period; throughout the present work, we fixed r = 4.0 
and A = 8.0 in units of a by simulating the geometry 
of actual peanut-shaped Cgo polymers whose geometry 
is reproduced by imposing a = 1 A. We found that U n 



1.5 



1.0 - 



0.5 



0.0 
-0.2 



5r/a = 1 .0 
Ala = 9.12 




n=3 



n=2 

n=1 
n=0 



10 



12 




4(z)/a 



£(z)/a 



FIG. 2: (color online) Profiles of the curvature-induced effec- 
tive potential Un{£) for one period [0, A]. Geometric param- 
eters ro = 4.0 and A = 8.0 in units of a are fixed. Integers 
n represent the angular momentum of eigenstates in the cir- 
cumferential direction of a hollow tube. 



takes extrema at £ — (or A) and £ = A/2, where r takes 
the maximum (r = ro) and the minimum (r = ro — Sr) 
values, respectively. It also follows that U n for n = is 
negative for any £ as derived from the definition of U n . 
To solve Eq. ([?]), we use the Fourier series expansions 



Un® = E 

G = ±2nj/A (j 
pansions into Eq 
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where 



= 0, 1, 2, • • • ). Substituting the ex- 
([4]), we obtain the secular equation 

Eg=-g c U G )c k-G = that holds for 
all possible ks and ns. The summation has been trun- 
cated by G c = 20-71-/ A because of the rapid decay of Ug 
with \G\. We then numerically calculated eigenvalues for 

< k < 7r/A and evaluated the low-energy band struc- 
ture for several different Sr, as depicted in Fig. [3l In all 
the cases, there is some energy gap at the Brillouin zone 
boundary Go = tt/A, where a wider energy gap occurs 
for a larger Sr, as expected from the large amplitude of 

1 [/„(£) I with increasing Sr (see Fig. [2]). 

We now consider the Coulombic interactions between 
spinless electrons. The interactions make the electron- 
hole pairs share the ground state of the noninteracting 
electron system, wherein the most strongly affected states 
are those lying in the vicinity of Ep. As a consequence, 
the single-particle density of states n(u>) near Ep exhibits 
a power-law singularity of the following form (34j 



n(to) on \hu> — Ep\ a , a 



K 



1. 



(5) 



The explicit form of K is derived using the bosonization 
procedure [34[ as follows 



Here. 

54(g) 



K=JhnJ ^* + "®-»® (6) 
9^0 y 2nhvp + g A (q) + g 2 (q) 

vp — h~ 1 dE/dk\k=k F is the Fermi velocity, and 
= V(q,m) and g 2 (q) — V(q,m) — V(2kp,m) 
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FIG. 3: (color online) Energy-band structure of a sinusoidal 
hollow tubule with ro = 4.0 and A = 8.0 in units of a. 



are ^-dependent coupling constants. V(q,m) is the 
Fourier transform of the screened interaction V(r) = 
— e 2 e _K ' r '/(47re|r|), where eis the dielectric constant and 
k is the screening length. The transformation is per- 
formed in terms of the curvilinear coordinates (£, 9), and 
thus, the resulting V(q,m) becomes a function of both 
the momentum and angular momentum transfers, q and 
m, respectively. To make concise arguments, Ep was as- 
sumed to lie in the lowest energy band (n = 0). This 
allows us to eliminate the index m from V(q,m), which 
leads to V(q) = -j^log[(g 2 + K 2 )r%] for qr <C 1, in 
which ro serves as the short-length-scale cut-off. 

The aim of the present study is to examine the Sr- 
dependence of K and a, which requires quantification 
of kp and vf- Among the many alternatives, we set 
kpa = 2.0£ x 10 -2 (£ — 1, • • • , 8), and evaluate the dimen- 
sionless Fermi velocity defined by vp = de / d{ka)\k = k F 
for each kp] the original vf is recovered by the relation 
vf = hvp/(2m*a). We chose the maximum value of kp 
such that it does not to exceed the <5r-dependent zone 
boundary Go for all Sr; it readily follows that Go takes 
the minimum value at Sr = ro, resulting in Goa ~ 0.17, 
and hence, kp < Go can be satisfied for all Sr. 

Figure [U (a) shows the plot of de/dk as a function of 
ka, from which vp for various kp and Sr are deduced, as 
shown in Fig.|3](b). Figure0](b) shows that Vp for all kp 
remains constant for Sr/a < 2.0, but decreases remark- 
ably with increasing Sr for Sr/a > 2.5. The decrease in 
vf is caused by the decrease in the slope of dispersion 
curves e = e(k) with Sr (see Fig. [3]). The results shown 
in Fig. 2] (b) imply that a change in Sr leads to a quanti- 
tative alteration in K and a, particularly at Sr/a > 2.5, 
and furthermore, the degree of alteration increases for 
larger kp, as will be demonstrated later. 

Figure [5] shows the 5r-dependence of both K and a 
for different kp values. According to the bosonization 
procedure [34| , we set the screening parameter k to be 
Ka = 1.0 x 10 -3 , which is smaller than all kpa values 



FIG. 4: (color online) (a) fc-dependence of the derivative 
de/dk for the lowest energy band, (b) <5r-dependence of the 
dimensionless Fermi velocity vp = de/d(ka)\k=k F for different 
kpa with increments A(kpa) — 0.02. 



that we have chosen. We also set the interaction-energy 
scale e 2 /(47rea) to be 1.1 in units of h 2 /J2m* a 2 ) by sim- 
ulating that of Ceo-related materials [43, Ej . The in- 
sets in Fig. [5] shows the ^-dependence of K and a at 
Sr/a — 2.0; each K and a take the minimum and max- 
imum values, respectively, at the specific kp satisfying 
the relation (d/dkp) \og[(2Tr hvp + <?4)/<72] = 0, which is 
equivalent to dK/dkp = 0. This result is the same for all 
Sr at Sr/a < 2.0. 

The salient features of Fig. [5] are the significant de- 
crease in K and increase in a with an increase in Sr for 
Sr/a > 2.5, as predicted earlier. It is interesting to note 
that such 5r-driven shifts in K and a are attributed to 
the effects of geometric curvature on the nature of TLL 
states. In fact, an increase in Sr amplifies the curvature- 
induced effective potential U n (£), thus yielding a mono- 
tonic decrease in vp at Sr/a > 2.5 (see Fig. @1. The 
decrease in vp plays a dominant role in the numerator 
of the expression in Eq. ([6]), and eventually leads to the 
systematic shifts in K and a. We have confirmed that a 
change in the value of k does not substantially affect the 
behaviors of K and a in a qualitative sense, while the 
absolute values of K and a moderately depend on the 
choice of k. 

An important consequence of the results shown in 
Fig. O is that nonzero surface curvature yields diverse 
alterations in the TLL behaviors of deformed cylinders. 
This is because various kinds of power-law exponents 
observed in TLL states are related to the quantity K; 
some such exponents are the exponent of power-law de- 
cay in Friedel oscillation 41| and that of temperature (or 
voltage)-dependent conductance [42j]. The present theo- 
retical predictions need to be confirmed experimentally, 
thus opening a new field of science that deals with quan- 
tum electron systems on curved surfaces. 

In conclusion, we reveal that the power-law exponent 
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FIG. 5: (color online) <5r-dependences of K and a denned by 
Eqs. ((5]l and ©, respectively. The screening parameter k is 
set to be kcl — 1.0 x 10~ 3 . Insets: Nonmonotonic behaviors 
of K and a as a function of kp at Sr/a — 2.0. 



a of the TLL states in deformed hollow nanocylinders 
shows a monotonic increase with an increase in the de- 
gree of surface curvature. The increase in a is attributed 
to the curvature-driven effective potential f/ n (0 that 
acts on electrons moving along the curved surface. The 
present results suggest that there are shifts in the power- 
law exponents of TLL states of real low-dimensional 
materials such as the peanut-shaped Cgo polymers and 
MoS 2 . 
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